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Apparatus and method for the processing of sectional images 



5 The invention relates to a method and an apparatus for the processing of 

sectional images that are reconstructed from X-ray projections of an object, preferably 
from cone beam projections. 

Reconstructed 3D sectional images in C-arm based volume imaging 
often exhibit spatially slowly varying inhomogeneities caused by inconsistent 

10 projection data due to various reasons. First, scattered radiation constitutes one of the 
main problems. Especially for system geometries with large cone angle and therefore a 
large irradiated area, such as C-arm based volume imaging, scattered radiation produces 
a significant background added to the desired detected signal and thus significant 
cupping shaped inhomogeneities occur in the reconstructed volume. Further sources of 

1 5 inhomogeneities are beam hardening (cupping) and insufficient truncation correction 
(cupping and inverse cupping). 

When a tight gray level window is applied to the viewed slices in the 
reconstructed 3D image in order to find low contrast details, the inhomogeneities cause 
parts of the region of interest to be either below or above the edges of the applied window. 

20 Low contrast details are therefore much more difficult to find in presence of 
inhomogeneities. 

Anti-scatter grids reduce the amount of detected scattered radiation, but 
the reduction is not sufficient for typical system geometries for volume imaging, by 

25 additionally increasing the signal-to-noise ratio. Different algorithms for projection 
based scatter compensation have been proposed (cf. Maher K. P., Malone J. F., 
"Computerized scatter correction in diagnostic radiology", Contemporary Physics, vol. 
38, no.2, pp. 131-148, 1997), but accurate quantitative scatter estimation is difficult to 
achieve. Truncation correction schemes like elliptical extension of projections (as 

30 described in R. M. Lewitt, Processing of incomplete measurement data in computed 
tomography, Med. Phys., vol. 6, no. 5, pp. 412-417, 1979) can prevent severe shading 
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artifacts at the boundaries of the FOV, but do not completely correct for global shading. 
Beam hardening correction is possible, but requires a deep knowledge of acquisition 
parameters, that is not always available. 

Based on this situation it was an object of the present invention to 
5 provide simple means for improving sectional X-ray images that are affected by 
inhomogeneities. 

This object is achieved by an apparatus with the features of claim 1 and a 
method with the features of claim 10. Preferred embodiments are subject of the 
1 0 dependent claims . 

The apparatus according to the present invention is suited for the 
processing of sectional images of an object, wherein such sectional images are 
reconstructed from X-ray projections of the object that were taken from different 
directions. The apparatus may optionally be adapted to execute the reconstruction of the 
1 5 sectional images from projections or it may alternatively receive already reconstructed 
sectional images for further processing. The apparatus is adapted to execute the 
following steps with a sectional image under consideration: 

The determination of a baseline function that at least 
approximately describes spatially slowly varying artifacts of the 
20 image. As will be explained in more detail below, such kind of 

artifacts are quite typical for a class of sectional images 
reconstructed from X-ray projections. 

Calculating a corrected sectional image by compensating the 
original sectional image with the help of said baseline function, 
25 e.g. by subtracting the baseline function from the original image. 

The apparatus has the advantage to (approximately) correct certain 
slowly varying artifacts, wherein the correction particularly makes the inspection of a 
small gray- value window possible. Because the apparatus does not try to correct the 
artifacts on the level of individual projections but processes the whole reconstructed 
30 sectional image, it has the further advantage to operate quite fast. 

The sectional image that is processed by the apparatus may be a 
two-dimensional section through an object like the human body, or it may be three- 
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dimensional and represent a volume of interest in an object. In both cases, the baseline 
function is of the same dimensionality as the original sectional image. In case of a 
three-dimensional sectional image, the baseline function may be composed of separate 
two-dimensional baseline functions, wherein each of these two-dimensional baseline 
5 functions is calculated for a two-dimensional slice from the original three-dimensional 
sectional image. In this case the separate two-dimensional baseline functions are 
preferably fitted to each other in neighbouring slices in order to achieve a smooth 
function all over the whole volume. 

According to a preferred embodiment of the apparatus, the determination 
10 of the baseline function comprises the following two sub-steps: 

The segmentation of areas from the original sectional image, in 
which the reconstructed X-ray density lies within a given 
interval. This interval may especially comprise the X-ray density 
of water, such that the segmented areas correspond to regions of 
1 5 the object that primarily comprise water. 

Determination of the baseline function based only on data that 
correspond to said segmented areas. Preferably the segmentation 
is followed by the operation of erosion in order to avoid 
misleading results at the boundaries of the segmented areas. 
20 The restriction of the search for the baseline functions to certain areas of 

the sectional images is especially successful in images of biological objects like the 
human brain. In this case, the image can roughly be subdivided into three classes, i.e. 
bone, air, and soft tissue, the latter corresponding in its properties approximately to 
those of water. The tissue areas may then be taken as a reference region for the desired 
25 calculation of the baseline function. 

In a preferred embodiment of the apparatus described above, the baseline 
function is determined by fitting a parametric model function to the data in the 
segmented areas. The model function may be especially a polynomial (preferably of 4th 
to 6th degree) and/or a spline function. 
30 According to another embodiment of the apparatus that calculates 

segmented areas, the baseline function is determined by low-pass filtering of data in the 
segmented areas. This approach is especially suited in cases where the segmented areas 
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are spatially distributed over the whole image area such that the low-pass filtering will 
be able to fill the gaps in between them. 

Still another approach for the determination of the baseline function is 
based on spectral filtering. In this case, the apparatus is adapted to execute the 
5 following steps: 

A spectral analysis (e.g. the three-dimensional Fourier 
decomposition) of the sectional images, or, in the case of an 
apparatus that segments the image, of the segmented areas of the 
sectional image. 

10 A composition of the baseline function from only the lower 

frequency components of the resulting spectrum. 
An advantage of such an apparatus is that it may make use of existing 
algorithms and software libraries for spectral processing of images. 

According to another embodiment of the apparatus, image areas outside 
15 the object are segmented in the original sectional image, and no correction of the 

original image is done there (i.e. the baseline function is set to zero in these areas). By ■ 
this way stability of the results is improved. 

The apparatus may especially comprise a rotational cone beam X-ray 
device for the generation of X-ray projections of an object. These projections may then 
20 be used for the reconstruction of the sectional image that is further processed by the 
apparatus. Especially in cone beam X-ray projections, there is a considerable 
contribution of inhomogeneous artifacts that need to be corrected. 

The invention further comprises a method for the processing of a 
sectional image that is reconstructed from X-ray projections of an object from different 
25 directions, comprising the following steps: 

a) The determination of a baseline function that describes spatially 
slowly varying artifacts of the sectional image. 

b) The calculation of a corrected image by compensating the 
sectional image with the help of said baseline function. 

30 The method comprises in general form the steps that can be executed 

with an apparatus of the kind described above. Therefore, reference is made to the 
preceding description for more information on the details, advantages and 
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improvements of the method. 

In summary, the present invention proposes a retrospective 
homogenization procedure. In a preferred embodiment, the reconstructed volume is 
segmented into bone, tissue and air based on their gray values, and the tissue-regions 
5 are used as support in order to fit a spatially slowly varying baseline function 

representing the smooth shape of cupping or other inhomogeneities. Finally the inverse 
of the estimated baseline is subtracted from the original slice to correct for the 
inhomogeneities. 

Beside fast calculation another advantage of the proposed method is that, 
10 in contrast to projection based correction schemes, the retrospective homogenization 
does not necessarily need detailed knowledge of system or acquisition parameters. 

These and other aspects of the invention will be apparent from and 
elucidated with reference to the embodiment(s) described hereinafter. 

In the following the invention is described by way of example with the 
15 help of the accompanying drawings in which: 



Fig. 1 is a diagram of the homogenization algorithm according to the 
present invention; 

20 Fig. 2-5 depict the processing of a simple example slice in different 

stages of the homogenization algorithm. 



Figure 2 depicts a simple sectional image I through a cylindrical test 
25 volume that comprises deeply dark and light structures (representing air and bone in a 
human body, respectively) embedded in a "tissue" of a mediate gray value. Moreover, 
there are some structures of small contrast with respect to the "tissue" and a spatially 
slowly varying inhomogeneity. In the lower part of Figure 2, a diagram of the gray 
value distribution along horizontal line through the image (at y = 200) is shown. When 
30 a tight gray level window is applied to the slice of the reconstructed 3D image in order 
to find low contrast details, the inhomogeneities cause parts of the region of interest to 
be either below or above the edges of the applied window. Low contrast details are 
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therefore much more difficult to find in presence of inhomogeneities. 

The basic idea of the present invention is to first find in each slice I of a 
3D volume a spatially slowly varying 2D baseline representing the smooth shape of 
cupping or other inhomogeneities and second to subtract this baseline from the original 
5 image. 

In order to find the baseline, two main problems have to be solved. First, 
a set of voxels has to be determined that can be used as supporting points for the 
baseline, and second, a good fit through these points has to be found. An appropriate 
selection of supporting points can be found by segmenting all voxels into three main 

10 classes of matter based on their gray value: air, tissue, and bone. Since the human body 
consists mainly out of tissue with a linear attenuation coefficient close to water, this can 
be done by applying a threshold around the linear absorption coefficient of water. The 
value of the linear absorption coefficient of water can be obtained by calculating it 
based on the applied X-ray spectrum or by statistical methods applied to the 3D images. 

15 To find a fit through all irregularly distributed supporting points, several 

methods can be used e.g. based on low-pass filtering, polynomial, and spline fitting or 
fitting an arbitrary parametric function. The main requirements are to globally model a 
smooth shape and to not follow local edges. By the choice of the function, application 
specific prior knowledge can be taken into account. Such information could be, e.g., the 

20 typical shape, area or gray level intensity of the inhomogeneities, that can be derived 
from previous example fits specific for the application. Finally, the estimated cupping 
baseline is subtracted from the original slice to correct the inhomogeneities. One of the 
advantages of the proposed retrospective homogenization algorithm is that it does not 
necessarily need detailed knowledge of system or acquisition parameters. 

25 An implementation example of the proposed homogenization algorithm 

is sketched in Figure 1 . Images illustrating the different steps of the algorithm are 
shown in Figs. 2-5 for a simple mathematical example slice. 

The algorithm may be carried out by a computer 1 being connected to an 
imaging apparatus like a cone beam CT (not shown). The computer will be provided 

30 with software for executing the algorithm which is illustrated by the flow chart of 
Figure 1. 

In block 10, the algorithm receives as input the original sectional image I 
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(e.g. a slice of a three-dimensional image). By taking into account the cupping of the 
original reconstructed volume, a threshold window between -200 HU and +50 HU is 
applied to the image I in block 1 1 . Since at the edges between tissue and bone or tissue 
and air steep transitions of the gray values appear that are partly within the threshold 
5 window but are not beneficial for the baseline estimation, these border points are 

eliminated by a morphological erosion algorithm in block 12. Starting from the original 
slice I, the resulting mask M of voxels selected as "valid tissue" is thus achieved in 
block 13 and shown in Figure 3 for the example slice of Figure 2. 

A fit through all supporting points of mask M is calculated in block 14 

10 by fitting a two-dimensional polynomial of sixth order to the data. The obtained 

baseline B of block 15 is shown in Figure 4 for the example slice of Figure 2. Despite 
the high order of the polynomial function, stable results are achieved due to the large 
number of available supporting points. In order to avoid stability problems of the fit 
outside the imaged object, the region of the object is segmented in block 16 by applying 

15 a second threshold followed by a four-neighbor connectivity algorithm. Finally, the 
estimated baseline is subtracted from the original slice to correct for the 
inhomogeneities, and the gray level of the image may be shifted to an arbitrary target 
value like the gray value of water according to the formula 

P(x,y) = I(x,y)-B(x,y) + n 

20 with |j, being a predetermined constant, e.g. representing the gray value 

of "water". The result is the corrected image I* in block 17 or Figure 5 for the example 
slice of Figure 2. Note that the homogenization algorithm is not intended to restore 
exact Hounsfield units. Rather, its purpose is to improve the detectability of low 
contrast details to facilitate the application of tight gray level windows. 

25 A similar procedure can also be applied in 3D space. The only difference 

would be to fit the baseline function not for each 2D slice separately, but a joint fit 
directly in 3 dimensions. Another extension of the described method compares the 
results of neighboring slices demanding for similar functions for the shape of the 
background. 

30 The estimation of the smooth shape of cupping alternatively can also be 

done be means of a Fourier decomposition and taking into account only the few first 
Fourier components. Depending on the data this can be done without the described 
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segmentation or applied on the segmented image M. 

Finally it is pointed out that in the present application the term "comprising" does not 
exclude other elements or steps, that "a" or "an" does not exclude a plurality, and that a 
single processor or other unit may fulfill the functions of several means. Moreover, 
5 reference signs in the claims shall not be construed as limiting their scope. 



